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Two-particle two-hole (2p2h) effect on the Gamow-Teller (GT) transition for neutron-rich nuclei is 
studied by the second Tamm-Dancoff approximation (STDA) with the Skyrme interaction. Unstable 
and ^'^Si, and stable "^^Ca nuclei are chosen to study the quenching and fragmentation of the 
GT strengths. Correlation of the 2p2h configurations causes 20% quenching and downward shift 
of GT giant resonances (GTGRs). The residual interaction changing relative angular momentum, 
appeared in the tensor force part, gives a meaningful effect to the GT strength distributions. In this 
work, 17 to 26 % of the total GT strengths are brought to high energy region above GTGRs. In 
particular, the tensor force brings strengths to high energy more than 50 MeV. STDA calculation 
within a small model space for 2p2h configuration is also performed and experimental data of '*®Ca 
is reproduced reasonably. 
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I. INTRODUCTION 

The Gamow-Teller (GT) transition is a suitable probe 
to investigate nuclear spin-isospin responses as well as 
the information on spin dependent channels of the nu¬ 
clear force. It also contributes to matrix elements of the 
nuclear /3-decay and double /3-decay. The former has to 
do with the time scale of rapid neutron-capture process 
(r-process) [i| and the prediction of decay heat of fission 
products important for the reprocessing of nuclear waste. 
The latter gives an important insight on the mechanism 
of neutrino mass @|- Low-lying GT state also points out 
the existence of SU(4) “supermultiplet” symmetry in nu¬ 
clei d, Hi¬ 
lt is well known that the total sum of GT strengths 
follows the model independent Ikeda sum rule Q, which 
is defined as S- — S+ = 3{N — Z). Here, S± denote 
the total sum of GT strength for /3^ transitions. How¬ 
ever, experimentally observed GT strengths around GT 
giant resonances (GTGRs) are only 50 — 60% of the to¬ 
tal value for the wide range of nuclei dQ. The leading 
origin of this damping (so-called quenching) mainly lies 
in the coupling with A-hole state or the correlation with 
higher-order configurations (see Refs. , and ref¬ 

erences therein). A recent experiment implies that the 
mixing of second-order configuration such as two-particle 
two-hole (2p 2h) state is more important than A-hole cou¬ 
pling §, lll| . The mixing of 2p2h states also accounts for 
the spreading width of GT stre ngths, which is studied by 
various theoretical approaches (l2 - ll3 - 

In the higher-order configuration mixing, the tensor 
force comes to play an important role, because it becomes 
more effective when the momentum transfer between two 
particles is high. Bertsch and Hamamoto studied the ef¬ 
fect of the tensor force on the total sum of GT strength of 
®°Zr including 2p2h states perturbatively and found that 
it plays a comparable role to the central force HI- Or- 
landini et al. obtained the same result in doubly-magic 


N=Z nuclei by using an energy weighted sum-rule ap¬ 
proach [ 23 . It was also shown that tensor force plays a 
important role in matrix elements of crT_|_ operator at low 
energies (see Ref. [IJ, and references therein). Drozdz 
et al. also studied the effect of the tensor force by sec¬ 
ond random-phase-approximation (SRPA) (^ . which in¬ 
cludes 2p2h states in a similar way to the standard RPA. 
However, their result showed a moderate effect of the 
tensor force in contrast to the other works. This differ¬ 
ence might be because they used a weaker tensor force 
than others. However, any of the above works does not 
consider the self-consistency in their formalisms which 
is essential to satisfy the Ikeda sum-rule, so that the ef¬ 
fects of the tensor force may not be represented correctly. 
Here, self-consistency means that one uses a same inter¬ 
action both in the ground state and the residual two-body 
interaction. In addition, they investigated excitation en¬ 
ergy up to at most 50 MeV. Bai et al. reported with a 
self-consistent RPA [23 that the tensor force can bring 
the GT strengths above 50 MeV already in a one-particle 
one-hole (Iplh) level. 

The 2p2h configuration mixing also affects the /3-decay. 
It is studied for light nuclei neighboring the /3-stability 
line [mill. In case of nuclei far from the /3-stability 
line, several author discussed the ,5-decays with the fi¬ 
nite rank separable approximation (FRSA) [H and the 
particle vibration coupling (PVC) [I!? which consider 
the 2p2h configuration mixing effectively via coupling to 
phonon states on the top of quasiparticle RPA (QRPA). 
The tensor force induces further modification of the low- 
lying GT state, and it is already studied in Iplh level 
[mil and FRSA [H- 

In this paper, we refocus the tensor force effect on 
the GT transition with the self-consistent Second Tamm- 
Dancoff approximation (STDA). As target nuclei, unsta¬ 
ble and ^"'’Si are chosen for future study of the 2p2h 
effect on /3-decay. We also investigate ^®Ca for com¬ 
parison with experimental data. STDA is, in fact, of 
a form omitting the ground state correlation part of Sec- 
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ond RPA (SRPA). However, it is not so significant in case 
of the GT transition so that STDA is suitable for 

the present purpose, because we are able to carry out it 
within a less computer resource than SRPA. In addition, 
STDA as well as SRPA is able to include a large model 
space and respects the Pauli principle which PVC and 
FRSA cannot take into account properly. 

The content of this paper is the following. Section II 
provides formalism of STDA and discuss the model space 
to be used in this work. Sec. Ill discusses our result of 
GT strength distribution and the quenching. In Sec. IV, 
the conclusion of present work and future plan are given. 


II. THEORETICAL METHOD 

The numerical calculation of the self-consistent STDA 
and SRPA had been difficult to be carried out until re¬ 
cently, but is now developed together with the evolution 
of computer resources and numerical technique, and is 
used to study monopole, dipole, quadrupole, and oc- 
tupole transitions up to ^'^Zr nuclei by realistic inter¬ 
action IlM as well as effective interactions such as 
Skvrme [32l - l34| and Gogny forces [35|. The present sta¬ 
tus of the self-consistent SRPA is summarized in Ref. 
( 3 ^ . We therefore illustrate our formalism used in this 
work briefly in sec. Ill A[ and discuss model spaces for 
STDA to be considered in sec. EBl 


A. STDA formalism 

STDA in this work uses single particle levels obtained 
by the Skyrme-Hartree-Fock (SHF) method [13 assum¬ 
ing spherical nuclear shape. We define the ground state 
of a nucleus as |SHF). SHF is solved in the coordinate 
space with a box boundary condition, rbox with a step 
size Ar = 0.1 fm. We use SGH Skyrme effective interac¬ 
tion for the central and spin-orbit forces [13, and Tel for 
the tensor force [13 • Variations of energy density and the 
spin-orbit potential by adding tensor force can be found 
in several papers, for example Ref. [23ll40l|- 

The basic formalism of STDA is same as SRPA formu¬ 
lated in Refs. [U H^, |4l|, 113 . An excited state |A) with 
respect to the ground state |0) is described by 

\\-JM) = Q\.j^\0), Qa;Jm|0) = 0 (1) 

where the phonon creation operator, Q\ are defined as 

^X;JM / j ^mi ^mi ^ ^mnijJpJh^rnnijJpJh’ 

mi 

( 2 ) 

The indices m and n denote particle states, while i and j 
denote hole states. The operators and 

create Iplh and 2p2h states coupled to the angular mo¬ 
mentum J and its projection to z-axis M, respectively 


[3l| - The two particle states and two hole states in the 
operator O are coupled to the angular momentum Jp 
and Jh, respectively. Compared with the phonon cre¬ 
ation operator of SRPA (3l|, the backward amplitudes 
characterized by and are omitted. In¬ 

stead, the ground state of STDA |0) is exactly identical 
to I SHF). 

We define the numbers of Iplh and 2p2h configura¬ 
tions in a given model space as Ni and N 2 , respec¬ 
tively. Only AT 2 = ±1 configurations are involved and 
non-charge exchange (AT^ = 0 ) and double charge ex¬ 
change (AT 2 = ± 2 ) configurations are decoupled from 
the phonon creation operator to reduce the dimension of 
STDA. 

The coefficients X and X are determined by solving 
the STDA equation, 


( A Ai2 
\^A21 A22 




(3) 


The JVi X JVi submatrix A is the standard Iplh RPA 
matrix. The submatrices A 12 (Xi x N 2 matrix) and A 21 
{N 2 X Ni) describe the coupling between Iplh and 2 p 2 h 
states, and A 22 (A ^2 x N 2 ) describes the coupling between 
2p2h states. The analytical forms are given in Ref. 

We calculate the matrix elements in the submatrices self- 
consistently, namely, the same interaction as the ground 
state is used. The rearrangement term appearing in ma¬ 
trices A [13 are also included. Since the matrix in the left 
hand of Eq. m is symmetric Hermite one, we can obtain 
eigenvalues, E\, directly by diagonalization, avoiding the 
imaginary solution problem, emerged in case of SRPA 
[^ Taking advantage of sparse matrices A, Eq. 

([5]) is solved by an appropriate numerical method. In this 
work, the EEAST linear algebra solver [d^l is used. 

If one omits the coupling between the 2 p 2 h states, A 22 
becomes diagonal and reads 


[A22] 
= (e 


mnij,m'n'i'j' 
m “t“ Cpi Ci 


j ')^mm' ^nn' 


(4) 


where the e are the single particle energies and x(m, n) is 
the antisymmetrizer between m and n states. This pre¬ 
scription is what is called the diagonal approximation, 
which reduces a computational task considerably. If this 
approximation works well, it would be useful for a qual¬ 
itative discussion of GT quenching. We will discuss the 
validity of the diagonal approximation in Sec. IHI PI 
Transition matrices of an operator F is given by 


(A; JM\F\0) = ^ -f ^ 

mi mnij 

where fmi = {rn\F\i) and fmnij = {mn\F\ij). In case of 
the GT transition, where F = 
strength function B±{GT) is written as 


2 


B^{GT) = 




( 6 ) 
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The 2p2h amplitudes X do not contribute B±{GT) in 
case of the one body external field as the present case. 
Hereafter, “GT” denotes f5~ type transition. 


B. Model Space 

Before carrying out the calculation of GT transition, 
we sought for an appropriate model space. In our calcula¬ 
tion, the continuum states are discretized by introducing 
a box boundary condition. Cutoff energy of unperturbed 
Iplh states is fixed to E].^^ = 100 MeV because it is 
insensitive to results. We also define the cutoff energy of 
unperturbed 2p2h state as 

First of all, we assessed the boundary box size denoted 
by 'f'box and a quantum number, N = 2n + l, where n and 
I are the number of node and orbital angular momentum 
of the single particle wave functions, respectively. The 
test was performed for and "‘^Ca under the condition 
of the diagonal approximation denoted by STDA(D), by 
setting ^ = 100 MeV. The results are shown in Fig. 
[T] and [5] for and "‘^Ca, respectively. The horizontal 
line corresponds to excitation energy of daughter nuclei. 
We used oi the AME mass table [d^. Let us begin 
with the upper panel (a) of Fig. [TJ This is the result 
of TV = 8 with different rbox- The strength distributions 
from E = —5 to 20 MeV do not converge for any rbox- As 
we increase iV, we found that the strength distributions 
for different rbox come close to each other. We finally 
obtain a good convergence for iV = 13 as shown in the 
bottom panel (b). The similar result is obtained for ^®Ca. 
In the upper panel (a) of Fig. the lines for different 
fbox do not converge for N = 8, while they become al¬ 
most the same shape for V = 12 as shown in the bottom 
panel (b). It means that it is enough to include a large 
N even in a relatively smaller rbox in order to obtain a 
reasonable convergence, at least within the present con¬ 
dition. We also notice from Fig[T] andthat, in contrast 
to the energy region from —5 to 20 MeV, the strength 
distributions at high energies are insensitive to rbox and 
TV. 

Next, we sought for appropriate E^^^. The test calcu¬ 
lations is performed altering E^^^ from 80 to 120 MeV by 
fixing TV = 13 for and TV = 12 for ^®Ga. The result is 
shown in Fig. [31 The GT resonances up to 10 MeV show 
a weak dependence on E'^^^ . Since we use the zero-range 
interaction in the residual interaction, it is not obvious 
how large E 2 p 2 h should be taken into account. In other 
words, we have to introduce an appropriate cutoff energy 
as long as we use the zero-range interaction. However, it 
would not be so rough to choose a particular E"^^^ from 
80 to 120 MeV in order to discuss the average behavior 
of the GT quenching, considering this weak dependence. 
We again notice that the GT strength distributions in a 
high energy region are less sensitive to E^^^. 

Finally, we adopt rbox = 10 fm, E^^^ = 100 MeV, and 
TV = 13 for and TV = 12 for “^^Ca, in this work. The 
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FIG. 1: (Color online) GT strength distribution of cal¬ 
culated by STDA(D) with various box boundary conditions 
from rbox = 8 to 14 fm. The upper and lower panels show 
the result for (a) A = 8 and (b) N = 13, respectively. The 
strengths are smeared by the Lorentzian function with a width 
1 MeV. is set to be 100 MeV. 



FIG. 2: (Color online) Same as Fig. [T] but for ^®Ca. The 
upper and lower panels show the result for (a) A = 8 and (b) 
A = 12, respectively. 


same model space as ^®Ga is adopted for ^^Si. Under this 
condition, the dimensions A 1 -I-A 2 are about 9x 10'^, 1.4x 
10*5 and 2.5 x 10^ for ^40, 34si, and ^SGa, respectively. 
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FIG. 3: (Color online) GT strength distributions of and 
^®Ca calculated by STDA(D) with different cutoff energies of 
unperturbed 2p2h states denoted by from 80 to 120 

MeV. 


III. RESULTS 

A. 2p2h effects on GT distributions 

We first discuss The GT strength distribution as 
a function of excitation energy of daughter nucleus from 
— 15 to 30 MeV is plotted in Fig. |4l The GT strengths 
are smoothed by the Lorentzian function with a width 1 
MeV, which represents the coupling to more complicated 
states. The position of experimentally observed 1+ state 
(1.8 MeV) is indicated by the arrow. The panel (a) illus¬ 
trates the TDA result. The GTGR appears at about 17 
MeV both for SGII and SGII+Tel. The low-lying reso¬ 
nances can be seen around 9 MeV for SGII and they are 
disturbed for SGII+Tel. In case of the STDA shown in 
the panel (b), the strength distribution of STDA is sys¬ 
tematically lower than TDA by about 6 MeV for SGII and 
the GTGR appears at about II MeV. This shift can be 
seen more strongly for SGII+Tel and several GT peaks 
appear at negative energies. 

The shift of resonances to lower energy region by the 
coupling with 2p2h states has been commonly observed 
and discussed in the other SRPA calculations [sol - lsl . . 
We will discuss this problem later in sec. mm 

Let’s return to Fig. 01 The peak height of GTGR 
for TDA is about 12. It is quenched to 8 for STDA 
as expected from 2p2h configuration mixing. TDA does 
not show any peaks around 1.8 MeV where the experi¬ 
mentally observed 1+ exists. On the other hand, STDA 
produces several resonances around this energy by shift- 



FIG. 4: (Color online) GT strength distribution of as a 
function of excitation energy with respect to its daughter nu¬ 
cleus. The upper, middle, and bottom panels are the results 
for (a)TDA, (b)STDA, and (c)STDA(D), respectively. The 
solid and dashed lines indicate the results calculated with 
SGII and SGII+Tel parameter sets, respectively. The ar¬ 
row indicates the experimentally observed 1'*' state [i^ . The 
strengths are smeared by the Lorentzian function with a width 
1 MeV. 

ing the low-lying GT resonances appearing at about 9 
MeV for TDA. The height of the low-lying resonance is 
not changed as much as GTGRs. This result is natural 
because GT states at low energies have less states to be 
coupled with them than those at high energies. There¬ 
fore, the contribution of the Iplh configuration to this 
state is still dominant. 

As seen in Fig. 01 GT strengths are distributed to other 
states by the 2p2h conhguration mixing. To see it more 
clearly, we plot the discrete GT strength in logarithmic 
scale in Fig. [5] The upper panel (a) and bottom panel 
(b) are the result of SGII and SGII+Tel, respectively. 
B(GT) of STDA are widely distributed, while that of 
TDA gives only several peaks around this energy region. 
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FIG. 5: (Color online) B(GT) of for (a)SGII and 

(b) SGII+Tel. The solid and dashed spikes indicate the re¬ 
sults of STDA and TDA, respectively. 

Figure [5] shows the GT distribution of ^^Si as a func¬ 
tion of excitation energies of its daughter nucleus from 
— 15 to 30 MeV, calculated by (a) TDA, (b)STDA and 

(c) STDA(D). Two experimentally observed I'*' states are 
denoted by the arrows, one of which is identical to the 
ground state of the daughter nucleus. The effects of 2p2h 
correlation and the tensor force are qualitatively same as 

The GT resonances of STDA are lower than TDA by 
about 7 MeV for SGII and roughly produce several reso¬ 
nances at which the observed I'*' states exist. SGII+Tel 
again shifts the GT resonances downward more strongly 
in case of STDA and produces the negative resonances. 
Similar to the heights and widths of low-lying reso¬ 
nances are also insensitive to the 2p2h effect. 

Figure [3 shows the GT strength distribution of "^®Ca. 
We plot experimental data measured by Yako et al. 
as well. TDA roughly reproduces the position of exper¬ 
imental GTGR for SGII and SGII+Tel. The low-lying 
resonance at 2.5 MeV is also reproduced fairly if we use 
SGII. However, widths of the resonances are not repro¬ 
duced at all due to the lack of coupling with higher-order 
configurations. STDA produces a slightly wider width, 
however, the position of GTGR seems rather low and the 
negative resonances appear as well as and ^^Si. 



FIG. 6: (Color online) Same as Fig. |4l but for ®^Si. 


action. The characteristic of the tensor force is to couple 
one state formed by two particles with another state dif¬ 
ferent in relative orbital angular momentum L by 2 as 
well as 0. Extracting the relevant angular momentum 
parts of the interaction matrix elements of tensor forces 
ytensor ^ emerged in A of Eq. ([3]), we obtain 


B. Tensor force effect and Appearance of negative 
GT resonance 

We observed that the GTGRs are significantly moved 
to lower energies by the tensor force in case of STDA. 
As a consequence, we obtained the unphysical negative 
resonances. In particular, the shifts induced by the ten¬ 
sor force was stronger than that by the central force. 
This would be because the tensor force becomes effective 
by high momentum transfer between two particles being 
possible. 

In order to investigate the tensor force effect more 
clearly, we pay attention to the role of the residual inter- 


V 'tensor \ '' 


(- 1 ) 


J' 


J' > 3fi 
1 jfl' jv' 




( L 1 J' 
L' 2 


G A 
1/2 1/2 1 
jfi A J' 


ly L' 'I 

1/2 1/21 L 

j^l' j^’ J' ] 


(7) 


where /i and v can be both particle and hole states, and 
and are the total and orbital angular momentums 
of the state /r (see (S^l for derivation of Eq. Q). 

We also consider the residual interaction of the central 
force, where there is a component interacting only among 
L = 0 or 2 states. The relevant angular momentum parts 
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FIG. 7: (Color online) GT strength distribution of ^®Ca as a 
function of excitation energy with respect to its daughter nu¬ 
cleus. The upper, middle, and bottom panels are the results 
for (a)TDA, (b)STDA, and (c)STDA(D) respectively. The 
solid and dashed lines indicate the results calculated with 
SGII and SGII+Tel parameter sets, respectively. The ex¬ 
perimental data is taken from Ref. [^ . The strengths are 
smeared by the Lorentzian function with a width 2 MeV for 
comparison with the experimental data as adopted in Ref. 
fs^ . which is the order of the experimental resolution. 


of the interaction matrix elements Reenter 

E ii,\\YL\\i,'){L\\YL\\L') 

L=0,2 

( L] ( I, I,, L] (8) 

X 1/2 1/2 1 W 1/2 1/2 1 L 

i Jm Jm' 1 J i > 1 J 

Note that the pairs of orbital angular momentum coupled 
to L (or L') are different in the center and tensor terms. 
Therefore, we consider them separately. We impose a 
simple condition that the interaction matrix elements of 
Eqs. ([7]) and ([5]) are zero unless the following case is 



HF 

case A 

case B 

case C 

case D 

TDA 

8.7 

15.4 

15.6 

9.7 

8.8 

STDA 

8.7 

11.2 

9.4 

5.9 

-0.9 

STDA(D) 

8.7 

10.3 

7.9 

5.6 

-0.7 


TABLE I: Center of mass of the GT peaks of defined as 
Ecm (see Eq. ®) in unit of MeV. 


satisfied; 

• case A : L = 0 for Eq.® (tensor part is excluded) 

• case B : L = 0, 2 for Eq.® (tensor part is excluded) 

• case C : L = L' for Eq. ® (central part is excluded) 

• case D : L = L'& ^ L' for Eq.® (central part is 
excluded). 

Case A retains the relative angular momentum of pairs 
[p/r'] and [I'ly'], and case B allows to exchange it by 2. 
Case C also retains the relative angular momentum of 
pairs [fiiy] and [pV'], and case D allows it to change. 

With the above conditions, we calculate the GT 
strength distribution for The center of mass of the 
GT peaks, defined as 


E 


cm 


<30MeV 

SieBiOOMeV 


(9) 


are also calculated and listed in Tab. m Figured] shows 
the results of the central part of the residual interaction, 
comparing with the Hartree-Fock (HE) result. For TDA 
shown in the panel (a), the GT resonances for case A 
appear at higher energies with respect to the HE result 
and form two sharp peaks. On the other hand, case B 
gives a minor change and shift Ecm upward only by 0.2 
MeV as seen in Tab. HI Namely, the GT resonances in the 
Iplh level is mainly produced for case A and the L = 2 
terms of case B are not of importance. However, we 
obtain the different result for STDA shown in the panel 
(b). Gase A doesn’t push the GT resonances as highly 
as TDA and case B shifts them to lower energies. 

Figure dj shows the results of the tensor part of the 
residual interaction. For TDA shown in the panel (a), 
the GT strength distributions of case C and D aren’t sig¬ 
nificantly different from the HE result. Ecm for case C 
and D are hence close to the one of the HE result, respec¬ 
tively, as shown in Tab. H] For STDA shown in the panel 
(b), the resonances of case C move to lower energies with 
respect to the HE result and Ecm is lowered by 2.8 MeV. 
Further downward shift can be seen for case D and Ecm 
is 6.8 MeV lower than the one of case C. In particular, the 
low-lying GT peak appearing around 1 MeV for case C 
shifts downward by as much as —10 MeV for case D and 
appears at around —9 MeV. The same result is obtained 
for other nuclei. Therefore, the residual interactions of 
the tensor force changing the relative orbital angular mo¬ 
mentum significantly contributes the downward shift of 
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FIG. 8: (Color online) GT strength distribution of cal¬ 
culated various conditions of the central part of the residual 
interaction for (a)TDA and (b)STDA(D). The thin solid line 
indicates the HF solution. The dashed and bold lines are 
cases A and B, respectively (see the text). 

the GTGRs. The center force part interacting between 
L = 2 states has also similar effect, however, its effect is 
not as strong as the tensor force. 

Besides the above discussions, one might think that 
the use of zero-range interaction contributes the strong 
shift of GT resonances in case of STD A. It is known that 
zero-range interactions such as the Skyrme force is not 
able to properly describe pairing correlation unless an 
appropriate cutoff energy is introduced [s^. Similarly, 
the use of zero-range interaction in STD A (SRPA) may 
induce a deviant value in 3 particle-1 hole or 1 particle-3 
hole matrix elements appeared in A of Eq. when 

high momentum transfers between particles (holes) oc¬ 
cur. If one uses a hnite-range interaction, a natural cutoff 
is able to be introduced. However, the downward shift is 
still commonly observed even when one uses finite-range 
interactions (^[ssLIS^I. 

It is recently pointed out that the downward shift 
obtained in the higher-order RPA is substantially at¬ 
tributed to the double counting between the residual 
interaction and the static response of the ground state 
(S^ , [ 5 ^. To overcome this problem, a subtraction 
method is suggested [12 • Gambacurta et al. applied 
it to SRPA and calculated monopole and quadrupole re¬ 
sponses. They showed that too low distributing strength 
functions obtained by SRPA are pushed up to as high en¬ 
ergy as RPA. It is therefore considered that the GT reso¬ 
nances appearing at low energy regions might be also due 



FIG. 9: (Color online) Same as Fig. |8l but for the tensor part 
of the residual interaction. The thin solid line indicates the 
HF solution. The dashed and bold lines are case C and D, 
respectively. 


to the same reason and they are expected to be shifted 
upward by adopting the subtraction method. 


C. Fragmentation and Quenching of GT strength 

We observed the fragmented GTGR in STDA and 
STDA(D) in the previous sections. Some of them are 
expected to be brought to higher excitation energies. To 
see them, we plot the GT strength distributions of 
^^Si, and "^®Ca in high energy regions above 20 MeV. 
Because of computational limits, we could not perform 
STDA at these energy regions so that only STDA(D) was 
carried out. 

Let’s start by looking at the result of shown in Fig. 
m For TDA shown in the panel (a), SGII does not give 
any significant peaks above 50 MeV, while SGII+Tel pro¬ 
duces a number of peaks there. This result is consistent 
with that pointed out in Ref. (^ . in which it is shown 
that the tensor force brings GT strengths to higher en¬ 
ergies even with Iplh RPA. In case of STDA(D) shown 
in the panel (b), we can see two marked results. The 
one is that SGII and SGII+Tel bring more strengths 
above 20 MeV than TDA. The second is that SGII+Tel 
significantly enhances the strengths in the energy region 
above about 60 MeV, while it suppresses them around 
30-45 MeV than SGII. Figures [TT] and [T^ plot the re¬ 
sults of ^^Si and "^®Ca at high energy regions, respec- 







































FIG. 10: (Color online) GT strength distribution of as 
a function of excitation energy with respect to its daughter 
nucleus. In contrast to Fig. [4l the horizontal axis is extended 
up to 100 MeV. The upper panel and lower panels are (a)TDA 
and (b)STDA(D) results, respectively. The resonances are 
smoothed by a Lorentzian function with 1 MeV width. 



FIG. 11: (Color online) Same as Fig. 1101 but for ^'^Si in the 
energy region from 20 to 85 MeV. 


lively. For both nuclei, STDA(D) brings a number of 
strengths to higher energy regions and SGII+Tel further 
enhances them above approximately 60 MeV and sup¬ 
presses around 30-40 MeV as well as 

To see the enhancement at high energies quantitatively, 
the sums of the GT strength of ^“^Si and "^^Ga in per¬ 
centages of the Ikeda sum rule are listed in Tabs. El uni 
and El respectively. We divide energy region into three 
parts, which are E < 20, 20 < E < 50, and 50 < E 
MeV. We also list the energy weighted sum rule, mi, for 
numerical check of the present code. Let’s first see the 
result of For SGII of TDA, 98.7% of the Ikeda sum 
rule is already exhausted in i? < 20. For SGII+Tel of 
TDA, it is reduced to 95.2%. For SGII of STDA(D), 



FIG. 12: (Color online) Same as Fig. 1101 but for ^®Ca in the 
energy region from 20 to 80 MeV. 


the sum of the GT strengths in A < 20 MeV is mean¬ 
ingfully quenched by 2p2h configuration mixing, being 
79.9%. SGII+Tel further promotes the quenching and 
the sum becomes 74.3%. For STDA, the quenching be¬ 
comes moderate as compared to STDA(D) and we obtain 
86.7% for SGII and 77.8% for SGII+Tel of the Ikeda sum 
rule in £1 < 20 MeV. 

Looking at the high energy regions above 20 MeV, a 
significant influence of the 2p2h configuration mixing can 
be seen. For STDA(D), totally 20.4% and 26.2% of the 
Ikeda sum rule are brought to A > 20 MeV for SGII and 
SGII+Tel, respectively. In particular, 8.0% of the Ikeda 
sum rule is found in 50 MeV < E for SGII and 17.1% for 
SGII+Tel. 

The same results are seen for ^^Si as well. In Tab. 
mil STDA as well as STDA(D) invokes the quenching 
in the energy region E < 20 MeV. STDA gives a some¬ 
what weaker quenching than STDA(D). The tensor force 
further promotes the quenching and the sums in A < 20 
MeV are reduced by 6.2%. In energy region 20 MeV < E, 
17.0% and 24.9% of the Ikeda sum rule are found for SGII 
and SGII+Tel, respectively, and we can still find mean¬ 
ingful GT sums above 50 MeV. 

For “^^Ga, the quenching below 20 MeV and enhance¬ 
ment above 20 MeV are qualitatively same as and 
^'^Si. The preceding work, the SRPA calculation [ 11 , 
estimates that 76% and 71% of the Ikeda sum rule are 
found in about E <20 MeV with and without the tensor 
force, respectively. These values overestimate the exper¬ 
imental data of 58.1% [111. In the present calculation, 
we obtain 88.5% and 80.3% in A < 20 MeV for SGII and 
SGII+Tel of STDA, respectively, which are about 10% 
larger than the SRPA calculation. On the other hand, 
the tensor force effect is closer to the SRPA calculation 
[H.!^ rather than Bertsch and Hamamoto’s result M- 
The present result is also close to that obtained in the 
PVG calculation [^, in case of SGII. 

It is known that non-energy (toq) and energy weighted 
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to 

0 

E<20 20<E<50 50<£' Total 

mi 

TDA 

SGII 

98.7 

1.4 

0.0 

100 

103 


SGII+Tel 

95.2 

3.4 

1.7 

100 

149 

STDA(D) 

SGII 

79.9 

12.4 

8.0 

100 

104 


SGII+Tel 

74.3 

9.1 

17.1 

101 

152 

STDA 

SGII 

86.7 

- 

- 

- 

- 


SGII+Tel 

77.8 

- 

- 

- 

- 


TABLE II: Sum of GT strengths of in percentages of the 
Ikeda sum rule for different energy regions. Energy weighted 
sum rule, mi, are also listed. 


3"Si 

E<20 20<E<50 50<E Total 

mi 

TDA 

SGII 

99.8 

0.6 

0.0 

101 

129 


SGII+Tel 

96.2 

3.6 

0.7 

101 

159 

STDA(D) 

SGII 

83.5 

10.0 

7.0 

100 

130 


SGII+Tel 

77.3 

6.8 

16.5 

101 

157 

STDA 

SGII 

88.2 

- 

- 

- 

- 


SGII+Tel 

79.8 

- 

- 

- 

- 


TABLE III: Same as TabUD but for ^'‘Si. 


sum rules (mi}of SRPA are analytically identical to 
those of RPA (sOl- It also holds for TDA and STD A. 
To see it, total GT sum (Total) equivalent to toq and mi 
are listed in Tabs. El uni and El We can confirm that 
TDA and STDA(D) give close values for all the nuclei. 


D. Diagonal Approximation 

In this section, we see how the diagonal approximation 
given in Eq. (j4]) works. Returning to Figs. |4l[6l and[7l we 
can see the GT resonance positions of STDA(D) are lower 
than those of STDA for all the nuclei. Except it, three 
interesting findings can be observed from these figures. 
The one is that the difference of GT resonance position 
between STDA and STDA(D) becomes smaller as we go 
from which is the lightest nucleus in this work, to the 
heaviest nuclei, '*®Ca. The second is that the difference 
between TDA and STDA(D) is much smaller in case of 
SGII+Tel than SGII. It implies that the diagonal matrix 
elements of the tensor force part would be larger than the 
off-diagonal ones. We can conhrm that it is not a wrong 
anticipation by looking at Ecm listed in Tab. H] For case 
B, which considers only the central residual interaction, 
the difference of Ecm between STDA and STDA(D) is 


«Ga 

E<20 20<E<50 50<E Total 

mi 

TDA 

SGII 

99.0 

1.7 

0.0 

101 

250 


SGII+Tel 

95.8 

4.2 

0.1 

101 

308 

STDA(D) 

SGII 

83.3 

11.3 

6.2 

101 

249 


SGII+Tel 

77.0 

8.1 

15.7 

101 

303 

STDA 

SGII 

88.5 

- 

- 

- 

- 


SGII+Tel 

80.3 

- 

- 

- 

- 


TABLE IV: Same as TabHH but for ^®Ca. 


1.5 MeV. On the other hand, if one considers only the 
tensor residual interaction, that is, case D shows only 
0.2 MeV difference between STDA and STDA(D). The 
third finding is that the overall GT resonance shape is 
not significantly different between STDA and STDA(D), 
for example, relative positions of low-lying GT resonance 
and GTGR. 

From the above results, we would say that the diago¬ 
nal approximation is a qualitatively good approximation 
when one calculates heavier nuclei or the tensor force is 
included. 


E. Smaller model space 

To take into account the 2p2h configurations as effec¬ 
tively as possible, a large model space has been intro¬ 
duced in the previous sections. However, the theoretical 
results do not necessarily reproduce experimental data, 
especially of "^®Ca. 

Toyama and Nakatsukasa studied giant dipole reso¬ 
nances in the TV = 82 isotones by the SRPA calcula¬ 
tion. They used a contact force in the residual interac¬ 
tion and a restricted model space smaller than the present 
study, and reproduced experimental data reasonably well 
[^ . PVC [ 2 I, m, m, [53 also reproduced experi¬ 
mental data successfully using a restricted model space 
of phonons coupling to Iplh states. An approach beyond 
RPA ( 2 ^ also adopts a relatively smaller model space and 
shows a good agreement with experimental data. 

We also performed a calculation with a smaller model 
space. As Toyama and Nakatsukasa did [sj, we con¬ 
sider the single particle levels near the Fermi energy of 
1 ^ 5 / 2 ) 1 '^ 3 / 2 j 2 si/ 2 j 1 / 7 / 212 p 3 / 2 ) 2 pi /2 and 1 / 5/2 orbits of 
proton and neutron to make the 2p2h states. The neu¬ 
tron and proton orbits are assumed to be fully occupied 
up to 1 / 7/2 a-iid 2 si/ 2 , respectively. The result is shown 
in Fig. [13l The downward shift of GT resonances is 
inhibited as compared with Fig. 0 GTGR appears at 
around 11 MeV which is close to the experimental one. 
The energy of the low-lying resonance is also reasonably 
reproduced. The obtained GT strength distribution is 
similar to that obtained in PVG [H^. The GT resonances 
of SGII+Tel locate at higher energy than SGII by about 
2 MeV. This is opposite to that obtained when the larger 
model space is considered. However, as we increase active 
single particle levels for 2p2h states, the GT resonances 
of SGII+Tel move to lower energies. The sum of GT 
strengths is also calculated and about 90% of the Ikeda 
sum-rule is exhausted for SGII and SGII+Tel in case of 
STDA, which is still higher than the experimental value. 


IV. CONCLUSION 

We studied the 2p2h effect on the GT strength distri¬ 
bution for ^"^Si and "“^^Ga. The strength distributions 
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FIG. 13: (Color online) Same as Fig. [7|in case of the restricted 
model space (see the text for the detail). 


were shifted to lower energy region by the 2p2h configura¬ 
tion mixing. The tensor force also induced further down¬ 
ward shift of the GT resonances. We showed that the 
tensor residual interactions changing relative orbital an¬ 
gular momentum contributes the strong downward shift. 
The central residual interaction exchanging relative or¬ 
bital angular momentum also played a similar role, but 
its effect was relatively small. 

One of the purposes of this paper was to see if we are 
able to describe the quenching of GT strengths by the 
self-consistent STDA. We have obtained the results that 
the sum of GT resonances below 20 MeV was quenched 
significantly in case of STDA. The tensor force further 
promoted the quenching of GT strengths and brought 
the missing ones to higher excitation energies. In spite 
of these results, however, we couldn’t reproduce the ex¬ 
perimental data of ^®Ga, and no improvements from the 
preceding SRPA (l8l.[^ was obtained. It was also found 
that the effect of the tensor force used in this work was 
similar to the result obtained in the preceding SRPA cal¬ 


culation, which gave a moderate tensor force effect than 
the perturbation approach fl^ . 

In the diagonal approximation, the GT resonances ap¬ 
peared at a lower energy region than full STDA system¬ 
atically and the resonances around GTGR were changed 
slightly. We found that the diagonal approximation rel¬ 
atively worked well, as we calculated heavier nuclei or 
included the tensor force. 

Following the method of Tohyama and Nakatsukasa 
[13, we calculated the GT distribution of "“^^Ga with the 
smaller model space. The obtained result reproduced 
the GT distribution of experimental data reasonably al¬ 
though the restricted model space of the single particle 
levels were chosen appropriately. At the moment, use 
of a rather restricted model space or cutoff of insignif¬ 
icant matrix elements as did in Ref. are the only 

approaches to apply STDA and SRPA to heavier or de¬ 
formed nuclei, in which much larger computer resource is 
usually required. The small model space would be there¬ 
fore helpful to obtain a physical insight of 2p2h effects as 
Tohyama and Nakatsukasa did [^ . 

Finally, we would like to note about the subtracted 
STDA and SRPA approaches (sec. IIIIBI) . Some of the 
GT resonances obtained in this work appeared at nega¬ 
tive energies, which are unphysical. However, the sub¬ 
tracted methods might have them shift upward and 
produce GT distributions at as high energy region as 
TDA(RPA). To check it is one of our future subjects. 
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